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TECHNICAL MEMORANDUM 


MULTIPLEXING READOUT CHANNELS IN PROPORTIONAL COUNTERS 


PROBLEM STATEMENT 

Proportional counters for x-ray astronomy under development at MSFC can use improved 
cathode planes constructed through a photolithographic process. Consequently, there is the 
opportunity to double the number of readout channels without incurring significant cost The 
signal from each channel must still be amplified before it can be processed. The problem is that 
there is no convenient way to make a corresponding increase in the number of amplifiers used in 
the instrument. This report describes and illustrates an idea for multiplexing pairs of channels 
into amplifiers, which allows the determination, with high degree of confidence, of the channel 
numbers and amplitudes that caused the signal. The instrument under consideration will have 
140 readout channels and 70 amplifiers in each of the x and y directions. 


DIFFICULTIES WITH MULTIPLEXING SCHEMES 

No multiplexing scheme can be completely successful. For example, suppose that an x-ray event 
interacting with the gas in the detector causes two distinct clusters of channels to produce signals, 
where one cluster consists of four adjacent channels (labeled using X’s below) and the other 
consists of two adjacent channels (labeled using Y). Each wire must be multiplexed with another 
wire in the array. The following situation can therefore occur. One of the wires in the 
four-channel cluster is multiplexed as shown: 

XXXX A 


The channel after A, B, must be multiplexed somewhere. If it is multiplexed with one of the 
wires from the two-channel cluster, 

XXXX AB Y 


it may be impossible to determine whether a photoelectron caused a two channel cluster AB or 
the legitimate two-channel cluster containing Y. Any multiplexing scheme is subject to this 
problem involving the creation of spurious two-channel clusters. 

Besides the occurrence of spurious two-channel clusters, multiplexing schemes are subject to two 
additional problems. A legitimate cluster can be enlarged when a channel adjacent to an 


endpoint of the cluster is multiplexed with a channel that is part of another cluster. Also, the 
amplitude of a channel that is part of a cluster can be increased if the channel is multiplexed with 
a channel that is part of another cluster. Both problems lead to error in the determination of the 
position of the event. In summary, any multiplexing scheme must face three problems: 

1. Spurious two-channel clusters 

2. Cluster enlargement 

3. Amplitude increase 


ASSUMPTIONS AND REQUIREMENTS 

Before describing the proposed multiplexing scheme, assumptions and requirements are given, as 
well as the general procedure for "decoding” the amplifier outputs. It is assumed that a 
legitimate x-ray event induces signals in two clusters of channels, where each cluster has width 
between two and five (although some clusters may be larger than five channels). We assume that 
if there were a continuum of channels, a cluster would appear as a Gaussian curve. In the 
detector within a single time slice (microseconds), one, two, or (rarely) three clusters of channels 
may produce signals. It is important to identify those situations where two clusters are activated. 

The general idea in reconstructing the actual channel clusters and amplitudes from the amplifier 
outputs involves first determining all the possible channels that might have been activated. Any 
channels that are "isolated" can be eliminated from further consideration ("isolated" means that 
neither adjacent channel was possibly activated by the event). For this idea to work, it is 
necessary that adjacent channels NOT be multiplexed with adjacent channels. Furthermore, 
channels that are multiplexed together must be separated by at least five channels. Otherwise, 
even a single cluster cannot be unambiguously reconstructed, as the following example shows: 

AXXXXB 


If A and B are both possible active channels that are multiplexed together, it is impossible to 
determine whether a five-channel cluster starts at A or ends at B. 


NUMERICALLY BASED MULTIPLEXING SCHEMES 

There are many approaches to obtaining a multiplexing scheme that will be good enough to allow 
the reconstruction of at least single clusters. The simplest idea is to multiplex channels in the 
first half of the array with channels in the second half of the array according to some numerical 
pattern. In this way, it is easy to guarantee separation of multiplexed pairs. With 140 channels in 
the array numbered from 0 to 139, we can consider the following scheme: 
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Given channel I for I between 0 and 69, multiplex channel I 
with channel k*I (mod 70) + 70. 

If k is chosen to be any integer relatively prime to 70, then the set of multiples of k (mod 70 ) 
will cover all of the integers between 0 and 69 (since k is invertible in the integers mod 70 if and 
only if k is relatively prime to 70). Note also that in this multiplexing scheme with k relatively 
prime to 70, given channel I for I between 70 and 139, the channel multiplexed with I is given by 
I*(l/k) (mod 70). For example, if k = 9, 1/k = 39 (mod 70), so the channel multiplexed with 
channel 88 is 88*39 = 2 (mod 70). Good choices for k are discussed below. 

The quality of the choice of k can be assessed using two measures. The first measure is the 
maximum number of consecutive channels in half of the array, which, if activated, can still be 
unambiguously reconstructed. This maximum unambiguously differentiable length (MUDL) is 
the number of "consecutive multiples" of k or 1/k (mod 70) that can be written in the integers 
mod 70 without obtaining any consecutive integers. For example, if k = 9, the progression of 
multiples 

0,9, 1 8,27,36,45,54,63,2, 1 1 ,20,... 

does not possess any consecutive integers until the 39th multiple, 1, is obtained. Notice, 
however, that for the multiples of 1/9 = 39 (mod 70), the progression proceeds 

0,39,8,47,16,55,24,63,32,1,... 

Thus for k = 9, in channels 0 through 69 MUDL is 39, while in channels 70 through 139 MUDL 
is 9. The MUDL is not simply k or 1/k (mod 70). For example, if k = 41, then k = 1/k (mod 
70) and it turns out that MUDL is 29. 

The other measure of quality for k involves determining the maximum number of consecutive 
channels spanning the center of the array, which, if activated, can still be unambiguously 
reconstructed. This central unambiguously differentiable length (CUDL) is simply 

min {b-a: (a,b) is a multiplexed pair and a < b}. 

For example, when k = 9 CUDL is 14 because (63,77) is a multiplexed pair. For k = 41, CUDL 
is 10 since 65 and 75 are multiplexed together. Use of the value computed for CUDL should 
perhaps be tempered by the determination of how many multiplexed pairs have the same CUDL 
and the next highest minimum separation. For example, though CUDL is 10 for k = 41, only one 
multiplexed pair is separated by 10. The next highest separation is 20. In contrast, for k = 29, 
CUDL is 14 with three pairs (58,72), (63,77), and (68,82) having the same separation, and the 
next highest separation is 28. 

In the schemes with k chosen to be 9 or 29 or 41, it is easily verified that any single cluster of 
width less than 9 can be unambiguously reconstructed from amplifier outputs. The best choice of 
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k to deal with pairs of clusters depends on characteristics of the detector’s operation. For 
example, if most clouds land near the middle of the detector, it would be best to use k = 9 and 
arrange to have channels 0 through 69 centered around the middle of the detector (since MUDL 
in that region is 39). If instead, clouds land in all portions of the detector with equal likelihood, it 
would be better to choose k to be either 29 or 41, because in each case k = 1/k (mod 70) and 
MUDL is 29. Then depending on the expected separation of clouds, a choice could be made 
based on CUDL. 


TECHNIQUES FOR RECONSTRUCTING ACTUAL CHANNELS 

With any of the numerically based methods described above, the reconstruction of actual 
channels and amplitudes involved in the event requires more than eliminating isolated channels. 
If we assume that triple events involving two-channel clusters are rare, it is possible to eliminate 
all occurrences of spurious two-channel clusters. The problem of cluster enlargement can be 
dealt with using a kind of smoothing. 


ELIMINATING SPURIOUS TWO-CHANNEL CLUSTERS 

The program called RECON.FOR in the appendix contains two subroutines used to eliminate 
spurious two-channel clusters. ELIM2A is invoked when a two-channel cluster is detected with 
two larger clusters. If both channels in the two-channel cluster are found to be multiplexed with 
channels that are part of the larger clusters, the two-channel cluster is eliminated. The reasoning 
is that triple events are rare to begin with, and the probability is remote that a triple event would 
occur with a two-channel cluster having both channels multiplexed into larger clusters. In fact, 
using actual data from a prior flight, we can obtain an empirical probability of 0.00518 for triple 
events. We can combine this with the probability that a channel will be multiplexed into one of 
two, four-channel clusters (8/138). Thus we obtain an estimate less than 0.00002 for the 
probability that a triple event has taken place and both channels of a two-channel cluster are 
multiplexed into the larger clusters. (This estimate should be reduced by the currently unknown 
conditional probability that one of the clusters of the triple event has width 2.) 

The subroutine ELIM2B is invoked when there are at least three clusters, at least two having 
width 2, and at least one having larger width. To handle this situation, we have the following 
result: 
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If a channel in a two-channel cluster is multiplexed under a 
numerical scheme as described above with k = 41 or k = 29 or 
k = 9, and its multiplexed pair is part of a cluster of width 3, 4, 
or 5, then the multiplexed pair of its neighboring channel 
cannot be adjacent to any of the other channels of the larger 
cluster or their multiplexed pairs. 

This situation is illustrated as follows: 

AAAAA CX 


The claim is that X cannot cause a spurious two-channel cluster with the multiplexed pairs of the 
largeT cluster. So, if there are spurious two-channel clusters, then any legitimate two-channel 
cluster would be multiplexed into only spurious clusters. Thus when invoked, ELIM2B 
eliminates all two-channel clusters containing a channel that is multiplexed into a larger cluster. 

The result described above can be proved by considering cases. One case is presented here; the 
others are similar. Suppose C > 69 is part of a two-channel cluster and is multiplexed (using k = 
41) with A, a channel that is part of a five-channel cluster. Then A = 41C (mod 70), and the 
channels multiplexed with C + 1 and C - 1 are (mod 70) 41(C + 1) and 41(C - 1), respectively. 
Clearly, neither of these can be equal to A + i or A - i for i = 1, 2, 3, or 4. Also, the channels 
multiplexed with C + 1 and C - 1 are less than 70, while channels multiplexed with A + i and A - 
i are greater than 70 (unless one cluster contains channel 70). Hence the multiplexed pairs of the 
channels in the two given clusters cannot be near each other to form a spurious two-channel 
cluster. (The case where one cluster contains channel 70 is verified by exhausting all the 
possibilities.) 


ALLEVIATING CLUSTER ENLARGEMENT 

To deal with the problem of cluster enlargement, subroutine ZAPEND identifies endpoints of 
clusters that have larger amplitudes than their neighbor channels. When a large amplitude 
endpoint is discovered, subroutine EXTRAPOLATE is called to determine what the expected 
endpoint value should be. This predicted value is computed by calculating the center of the 
cluster omitting the endpoint, looking at the amplitudes on the opposite end of the cluster, and 
using a linearly interpolated value. Exploiting the symmetry of the data in this way leads to a 
highly computationally efficient routine as well as excellent accuracy. The extrapolation 
subroutine was tested separately on a large set of clusters by comparing the actual cluster centers 
with the centers obtained using EXTRAPOLATE on those same clusters with an endpoint 
removed. The absolute difference between the predicted and actual centers of clusters had a 
mean of 4% of a channel width, and a standard deviation approximately equal to the mean. The 
maximum absolute difference was 46% of a channel width. A histogram showing the 
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distribution of the absolute differences is given in the appendix. Over 99% of the predictions are 
within 20% of a channel width of the actual center. Further testing results are given below. 


AMPLITUDE AMPLIFICATION 

The problem of amplitude amplification was not directly addressed. If the amplitude of an 
endpoint of a cluster turns out to be increased sufficiently, subroutine ZAPEND will adjust it 
appropriately. Errors in endpoint amplitudes have the largest significance in the computation of 
the event position. An amplitude increase caused by multiplexing in a central channel of a 
cluster obviously has negligible influence on the computation of the center of the cluster. But 
errors in multiplexed channels that are neither central nor on the end of a cluster lead to some 
error in the computed position of the center. Future work may involve assessing the magnitude 
of this type of error empirically, and perhaps designing a broader form of Gaussian smoothing to 
reduce the error. 


RESULTS 

The multiplexing scheme using k = 41 was tested using data derived from a previous balloon 
flight and randomly generated data. Single events were always correctly reconstructed, so they 
are not included in the results below. In addition, the program also correcdy reconstructed the 
event state (single, double, or triple event) in all cases. 

The detector in the previous balloon flight used 70 readout channels. The data from that flight 
were adjusted slightly to make it appear as if they came from a detector with 140 channels. This 
was done by doubling the beginning channel number of each cluster and building new clusters 
with the same widths as the old ones, but starting at the doubled locations. 

Testing was done by using existing MSFC analysis software to determine the centers of the 
events. Centers were first obtained for the actual data and placed in file CENTERSA. Then the 
actual data were supplied to a program that created file CODED consisting of the multiplexed 
amplifier outputs (file formats are given in the appendix). The reconstruction program, using 
only file CODED, built a file, PREDRAW, that contained the predicted channel numbers and 
their amplitudes. The analysis software was then run on file PREDRAW to obtain file 
CENTERSP. Finally summary statistics and a histogram were generated from the differences 
between CENTERSA and CENTERSP. One channel width on the 140 channel detector is 2 mm. 
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Using the 54 double events in the file from the previous balloon flight, 216 computations of 
cluster centers were made in each of CENTERSA and CENTERS P. The results are given below: 

Mean absolute difference = 0.024 mm 
Standard deviation = 0. 1 26 mm 
Maximum absolute difference = 1.28 mm 
Number of overestimates = 12 
Number of underestimates = 11 
Number of exactly (*) correct estimates = 193 

(*) to the limit of single precision arithmetic 

The file of simulated data contained 934 double events scattered uniformly on the detector. This 
gave rise to 3736 computations of cluster centers in each of the files CENTERSA and 
CENTERSP. The results are as follows: 

Mean absolute difference = 0.026 mm 
Standard deviation = 0.103 mm 
Maximum absolute difference = 1.28 mm 
Number of overestimates = 179 
Number of underestimates = 231 
Number of exactly (*) correct estimates = 3326 

A histogram shows these differences to be quite concentrated in the first 0.1 mm range: 


Bin 

Number of 

Range 

(mm) 

Observations 

0 - 0.1 

3456 

0.1 -0.2 

82 

0.2 - 0.3 

91 

0.3 - 0.4 

30 

0.4 - 0.5 

24 

0.5 - 0.6 

19 

0.6 - 0.7 

21 

0.7 - 0.8 

4 

0.8 - 0.9 

0 

0.9- 1.0 

5 

1.0- 

4 


Evidently multiplexing readout channels will not lead to significant error in the computation of 
x-ray event locations. 
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0.05 0.1 0.15 0.2 0.25 0.3 

Difference 

Figure A.l. Distribution of absolute differences between predicted and actual centers. 



APPENDIX B 


PRECEDING PAGE BLANK NOT FILMED 


File Formats 


RAW and PREDRAW 

Each event is associated with five records: 

(1) number of x channels number of y channels 

(2) x channel numbers active (range 1-140) 

(3) x amplitudes corresponding to active channels 

(4) y channel numbers active (range 1-140) 

(5) y amplitudes corresponding to active channels 


CODED 

Each event is associated with one record containing: 
number of x amplifiers active 
x amplifier numbers active 

x amplitudes corresponding to active amplifiers 
number of y amplifiers active 
y amplifier numbers active 

y amplitudes corresponding to active amplifiers 
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RECON. FOR 


(DECLARE 

PR06RAH RECON 
INTE6ER HAXCHAN 
PARAHETER (HAXCHAN=140 ) 

INTE6ER PX ,PY , XOUT , VOUT , «X , NY , PSUBX , PSUBY ,NXCLUHPS ,NYCLUHPS 
I NTE6ER XSTART , YSTART , XN I DTH , YN I DTH , XSUH , YSUH , I , J , N 
D I HENS I ON PX ( 0 : 69 ) ,PY ( 0 : 69 ) , XOUT ( 0 : HAXCHAN-1 ) , YOUT ( 0 : HAXCHAN-1 ) 
DIHENS I ON PSUBX (50), PSUBY (50) 

C 

C Reconstruct events in rax fora fro* events passing through the 
C preaaplifiers and coded in file CODED. The prograa proceeds by 

C assigning the observed amplitude froa a preaaplifier to both of 

C the channels aultiplexed to it. Then all isolated channels are 
C eliainated. An atteapt is aade to deteraine whether a single 
C event took place. If this cannot be established, subroutine 

C FIND2 is invoked, which is structured like an expert systea. 

C It contains rules which deal with duap enlargeaent and spurious 
C length 2 cluaps. 

C 

OPEN (J2, FILE* 'CODED' ,STATUS='OLD’ ) 

OPEN ( 33 ,FILE S ' PREDRAN ' , STATUS* 'UNKNONN' ) 

C Obtain preaaplifier levels for each record in CODED file 
C Beginning of aain loop 

C First initialize all preaaps to 0 since not all are in CODED 
10 DO 20 1*0,69 
PX(I)*0 
PY(I)*0 
20 CONTINUE 

DO 30 I*0,HAXCHAN-1 
XOUT ( I )*0 
YOUT ( I ) *0 
30 CONTINUE 
C 

READ(32,t,END*999)NX,(PSUBX(I),I*l,NX),(PX(PSUBX(I) ),I*1,NX), 

♦ NY, (PSUBY! I ) , 1*1 ,NY) , (PY(PSUBY( I ) ) , 1*1 ,NY) 

C 

C Convert preaaplifier outputs to all channels theoretically possible 
DO 40 1=0,69 
IF (PX(I).6T.O) THEN 
XOUT ( I )*PX ( I ) 

N=H0D(41$I, 701+70 
X0UT(N )*PX ( I ) 

END IF 

IF (PY(I).GT.O) THEN 
YOUT ( I )=PY ( I ) 

N=H0D(41II, 701+70 
YOUT (N1=PY(I ) 

ENDIF 
40 CONTINUE 
C 
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C Eliminate island* 

IF ( ( XOUT (0) . 6T.0) .AMD. ( XOUT ( i ) .LE.O) ) XOUT (O)--XOUT(O) 

IF (( Y0UT(0).6T.0I. AND. (YOUT(l). LE.O)) Y0UT(0)=-Y0UT(0) 

IF ( ( XOUT (HAXCKAN-1 ) . 6T . 0 ) . AND . ( XOUT ( HAXCHAM-2 ) . LE. 0 ) ) 

♦ XOUT ( HAXCHAM-1 ) =-XOUT ( ) 

IF ( ( YOUT (HAXCHAM-1 ) . 6T . 0 ) . AND . ( YOUT ( HAXCHAM-2 ) . LE . 0 ) ) 

♦ YOUT ( HAXCHAM-1 ) =-YOUT ( HAXCHAN-1 ) 

DO 50 1=1, HAXCHAM-2 

IF ((XOUT(I).6T.O). AND. (XOUT(I-l). LE.O). AND. 

♦ (XOUT(I*l).LE.O)) XOUT( I )=-XOUT( I ) 

IF ((YOUT(I). 6T.0). AND. (YOUT(I-l). LE.O). AMD. 

+ (Y0UT(I+1) .LE.O) ) YOUT ( I ) =- YOUT ( I ) 

50 CONTINUE 
C 

C Find the nuaber of active channel cluaps in each direction 
CALL COUNTCLUHPS ( XOUT , NXCLUHPS ) 

CALL COUNTCLUHPS l YOUT , NYCLUHPS) 

C 

IF (( NXCLUHPS. 6T.1). OR. (NYCLUHPS.GT.l)) THEN 
CALL FIND2I XOUT , YOUT .NXCLUHPS , NYCLUHPS ) 

ELSE 

C In this case, the nuaber of cluaps in each direction is 1 

C CLUHPINFO returns start channel nuabers, widths, and aaplitude 

C suas for the duap found 

CALL CLUf1PINF0( XOUT, 0, XSTART . XWIDTH, XSUH) 

CALL CLUHPINFO ( YOUT , 0 , YSTART , YN I DTH , YSUH ) 

IF ( (XMIDTH.NE. YHIDTH) .OR. (XSUH.NE. Y5UK) .OR. 

+ (XNIDTH.6T.5) ) THEN 

CALL F I ND2 ( XOUT , YOU T , NXCLUHPS .NYCLUHPS ) 

ELSE 

NRITE(33,' (213)' IXNIDTH, YHIDTH 
C Mrite out x channels triggered in raw foraat 
NRITE(33,81)(XSTART+J,J=1,XNIDTH) 

C Hrite out x aaplitudes 

«RITE(33,81 ) ( XOUT( XSTART+J ) , J=0, XtlIDTH-1 ) 

C Hrite out y channels triggered in raw foraat 
NRITE(33,81)(YSTART+J,J=1, YHIDTH) 

C Hrite out y aaplitudes 

HRI TE( 33 ,81 ) ( YOUT ( YSTART+J ) , J=0, YNIDTH-1 ) 

ENDIF 

ENDIF 

C End of aain loop 
60 TO 10 

81 FORHAT ( 10014) 

999 STOP 
END 
C 
C 

SUBROUTINE FIND2 ( XCHAN , YCHAN , NXCLUHPS , NYCLUHPS I 
C Top level processor for case where double events are suspected 
INTEGER HAXCHAN 
PARAHETER( HAXCHAN=140 ) 

I NTEGER XCHAN , YCHAN , NXCLUHPS , NYCLUHPS , I , INIT 

INTEGER HX(10) ,HY( 10 ) , XSTART( 10 ) , YSTART( 10) , XSUH( 10) , YSUH( 10) 

L06ICAL ELIH2A.ELIH2B 

DIHENS I ON XCHAN ( 0 : HAXCHAN- 1 ) , YCHAN ( 0 : HAXCHAN-1 ) 

DO 1 1=1,10 
XSUH(I)=0 
YSUH(I)=0 
1 CONTINUE 
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C Accumulate inforaation about the dumps 
INIT=0 

DO 10 1=1, NXCLUHPS 

CALL CLUHPINFO( XCHAN, INIT, XSTART (I),HX(I), ISUH( I ) ) 
INIT=XSTART(mNX(l) 

10 CONTINUE 
IN1T=0 

DO 20 I<1,NYCL(HB>S 

CALL CLUHP1NF01 YCHAN, INIT , YSTART II) ,NY U ) , YSUNl 1 ) ) 
INIT=YSTART(I)+NY(I) 

20 CONTINUE 

C Eliminate spurious endpoints of duaps with width 6T 2 
CALL 2 APENDI NXCLUHPS , XCHAN , XSTART , NX , XSUH ) 

CALL 2 APEND ( NYCLUNPS , YCHAN , YSTART , NY , YSUH ) 

C Sort the cluap info pointers by duap width 
30 IF (NXCLUHPS.6T.1) CALL SORT (NXCLUHPS, NX, XSTART, XSUH) 

IF (NYCLUHPS.6T.1) CALL SORT (NYCLUHPS, NY, YSTART, YSUH) 

C The follouing code operates like an expert system; each condition 
C of the case structure represents a known situation, and a Ion 
C confidence catch-all condition occurs at the end 

IF ( (NXCLUHPS. EO. 2). AND. (NYCLUHPS. EO. 2). AND. (NX(1).LT.HX(2)) 

4 .AND. (HX ( 1 ) ,EQ.NY( 1) ) .AND. ( HX ( 2 ) . EB . MY ( 2 ) ) ) THEN 
C Case uhere there are clearly 2 events of different Nidths 

CALL NR1TER(XCHAN, YCHAN, NXCLUHPS, NYCLUHPS, XSTART, YSTART, NX, NY) 
ELSE I F ( (NXCLUHPS . EQ . 2 ) . AND . ( NYCLUHPS . EB . 1 ) . ANO . 

♦ (NX(1).E0.NX(2)).AND. (NX(1) .EB.NY(l)) ) THEN 

C Case uhere events line up in the y dimension, but widths equal 

CALL NRITER( XCHAN, YCHAN .NXCLUHPS . NYCLUHPS , XSTART , YSTART , NX , NY ) 
ELSE IF ( ( NXCLUHPS . EB . 1 ) . AND . ( NYCLUHPS . EB . 2 ) . AND . 

♦ (HY(1).EQ.NY(2)) .AND. (NX(1) .EO.NY(l)) ) THEN 

C Case uhere events line up in the x dimension, but uidths equal 

CALL NRI TER ( XCHAN , YCHAN , NXCLUHPS , NYCLUHPS , XSTART , YSTART , HX ,NY ) 
ELSEIF ( ( NXCLUHPS . 6T . 2 ). AND .( NX ( NXCLUHPS- 1 ) . 6T . 2 ) . AND . 

♦ (HX(1).EB.2).AND. 

♦ (ELIH2A1 XCHAN, NXCLUHPS, NX, XSTART, XSUH)) ) THEN 

C Case uhere a uidth 2 duap is multiplexed with other duaps of 
C larger uidth. ELIH2A checks the duals, and if successful, reduces 
C the number of clumps, resets the uidth, start, and sum pointers, 

C and returns .TRUE. 

60 TO 30 

ELSEIF ( (NYCLUHPS. 6T. 2). AND. ( NY { NYCLUHPS- 1).6T. 2). AND. 

♦ (NY(1). EO. 2). AND. 

♦ (ELIH2AIYCHAN, NYCLUHPS, NY, YSTART, YSUH)) ) THEN 
60 TO 30 

C 

ELSEIF ( (NXCLUHPS. 6T. 2). AND. (NX (2) . EQ. 2) .AND. 

♦ (NX ( NXCLUHPS ).6T. 2) .AND. 

♦ ( EL I H2B ( XCHAN, NXCLUHPS , HX , XSTART , XSUH ) ) ) THEN 

C Case uhere a uidth 2 dump is multiplexed uith other uidth 2 
C dumps. ELIH2B looks for uidth 2 clumps multiplexed uith a larger 
C cluap and eliminates them. If any clumps are eliminated, ELIH2B 
C returns .TRUE. 

60 TO 30 

ELSEIF ( ( NYCLUHPS. 6T. 2). AND. (HY(2 ) .EB. 2). AND. 

♦ (ELIH2BIYCHAN, NYCLUHPS, NY, YSTART, YSUH)) ) THEN 
60 TO 30 

ELSE 

CALL NRITERIXCHAN, YCHAN, NXCLUHPS, NYCLUHPS, XSTART, YSTART, NX, NY) 
END IF 
RETURN 
END 
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LOGICAL FUNCTION ELIH2A(CHAN, NCLUHPS, NIDTH, START, SUN) 

INTEGER HAXCHAN 

PARAHETER(HAXCHAN=140) 

1NTE6ER CHAN, NIDTH, START, SUN, NCLUHPS, DUAL, l 
DIMENSION CHAN(0:HAXCHAN-1),NIDTH(10| ,START(10|,SUN(10) 

C Nhen there are aore than 2 cluaps in one diaension, deteraine 

C whether one of the width 2 cluaps is spurious by seeing if 

C the duals of its aeebers are part of other cluaps 

C 

ELIH2A=.FALSE. 

IF ( (MIDTHf 1) .EQ.2) .AND. 

♦ ( CHAN ( DU AL ( ST ART ( 1 ) ) ) . EQ . CHAN ( ST ART ( 1 ) ) ) . AND . 

♦ ( CHAN ( DUAL ( START! 1 ) +1 ) ) . EQ . CHAN (START (11+1)1 (THEN 
CHAN! START! 1 ) ) =-CHAN( START |1)) 

CHAN (START ( 1 ) +1 )=-CHAN ( START (1)+1) 

NCLUHPS=NCLUHPS-1 
DO 100 1=1 ,NCLUHPS 
WIDTH! I ) =MIDTH( 

START ( I )=START ( 1+1 ) 

SUM! I )=SUM( 1*1 ) 

100 CONTINUE 

ELIH2A=.TRUE. 

END IF 
RETURN 
END 
C 

LOGICAL FUNCTION EL IH2B ( CHAN , NCLUHPS , M I DTH , ST ART , SUN ) 
INTE6ER HAXCHAN 
PARAMETER ( HAXCHAN= 140 ) 

INTEGER CHAN , WIDTH, START ; SUN , NCLUNPS , DUAL , CLNO, I,J,K 
DIMENSION CHAN(0:HAXCHAN-1 ) , WIDTH! 10) , START ( 10) , SUN! 10) 
L06ICAL MENDER 

C If a width 2 duap is found which is aultiplexed with a larger 
C duap, it is eliainated. This function is only invoked when 
C there are at least 2 width 2 cluaps 
ELIH2B=. FALSE. 

CLN0=1 

10 IF ( (HIDTH(CLNO).EO. 2). AND. (CLNO.LT. NCLUHPS) ) THEN 
I =DUAL( START (CLNO)) 

J=DUAL ( START ( CLNO ) *1) 

C See if I or J is aultiplexed into the larger duap 
IF ( (NENBER(I, START! NCLUNPS ) .WIDTH! NCLUHPS I ) I . OR . 

♦ ( MEMBER ( J , START ( NCLUHPS ) , WIDTH! NCLUHPS ) ) ) ) THEN 

C Eliainate the current duap 

C Do not increaent CLNO in this case. 

CHAN! START ( CLNO ) ) =-CHAN (START ( CLNO ) I 
CHAN ( START | CLNO ) *1 ) =-CHAN ( STAR T ( CLNO ) ♦ 1 ) 

C Hove width, start, and sua pointers down 

DO 30 X=CLN0, NCLUHPS- 1 
HIDTH(IO=HIDTH(ICH) 

START! K)=START(K+1) 

SUH(K)=SUH(K+1) 

30 CONTINUE 

NCLUHPS=NCLUHPS-1 
ELIH28=.TRUE. 

ELSE 

CLN0-CLN0+1 
END IF 
60 TO 10 
ENDIF 
RETURN 
END 
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L06ICAL FUKCTIM HEHBER(CHNO, START, WIDTH) 

INTEGER CHNO, START, MIDTH 

C Determine whether CHNO is part of the duap of length WIDTH 
C beginning at channel nuaber START 

IF ( (CHNO. GE. START) .AND. (CHNO.LT.START+WIDTH) ) THEN 
RENDER 8 . TRUE. 

ELSE 

HEHBER 8 . FALSE. 

ENDIF 

RETURN 

END 

C 

SUBROUTINE NR I TER ( XCH AN , YCHAN . NXCLUHPS , NYCLUNPS , XSTART , YST ART , NX , NY ) 
INTE6ER HAXCHAN 
PARANETER ( HAXCHAN=140 ) 

INTEGER XCHAN , YCHAN , NICLUHPS , NYCLUNPS , TOTXCH , TOTYCH 
I NTE6ER NX ( 10) ,NY (10 ) , XSTART (10), YST ART ( 10 ) , I , J 
D I HENS I ON XCHAN(0 : RAXCHAN-1 ) , YCHAN ( 0: HAXCHAN-1 ) 

T0TXCH=0 

T0TYCH=0 

DO 10 I*1,NXCLUHPS 
TOTXCH=TOTXCH+NX(I) 

10 CONTINUE 

DO 20 1=1 (NYCLUNPS 
TOTYCH=TOTYCH+WY(l) 

20 CONTINUE 

NRITEI33,' (213)' 1T0TXCH, TOTYCH 
C Write out x channels triggered in raw foraat 

NRI TE ( 33 , SI ) ( (XSTART ( I ) + J , J*1 ,NX ( I ) ) , I =1 , NXCLUWS ) 

C Write out x aaplitudes 

NRITE(33,81)((XCHAN(XSTART(I)+J),J=0,WX(I|-1), 1=1, NXCLUHPS) 

C Nrite out y channels triggered in raw foraat 

NRITE(33,81)((YSTART(I)+J,J=1,WY(I)),I=1, NYCLUNPS) 

C Write out y aaplitudes 

NRITE(33, 81 )( (YCHAN(YST ART( 1 )+J),J*0,WY( I )-l),I*l, NYCLUNPS) 

RETURN 

81 FORHATI 10014) 

END 

C 

SUBROUTINE CLUHP INFO (CHAN, IN1T, START, WIDTH, SUN) 

INTE6ER HAXCHAN 
PARAMETER ( KAXCHAN-140 ) 

INTE6ER CHAN, INIT. START, WIDTH, SUN, I 
DIMENSION CHAN(OiRAXCHAN-l) 

C Find consecutive active channels and their aaplitude sua 
C starting search at position INIT 
C 

C Find the first non-zero aaplitude channel 
START*INIT 

5 IF (CHAN(START).LE.O) THEN 
START=START+1 

IF ( START. GE. HAXCHAN) THEN 
WIDTH*0 
SUH=0 
RETURN 
ELSE 
GO TO 5 
ENDIF 
ENDIF 
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C Find total width and accuaulate sua 
SUH=CHAN( START) 

NIDTH=1 

I=STflRT+l 

10 IF (CHAN(I).6T.O) THEN 
NIDTH=NIDTH+1 
SUH=SUH+CHAN(I) 

1=1+1 

IF (I.6T.HAXCHAN-1) RETURN 
60 TO 10 
END1F 
RETURN 
END 
C 

SUBROUTINE SORT ( NCLUHPS, N I DTH, START, SUN) 

INTE6ER NCLUHPS f HI DTH ( 10 ) , START ( 10 ) , SUN ( 10 ) , I , J 
C Sort the cluap width, start, and sua pointers by width 
C Bubble sort is appropriate here since NCLUHPS is 5 or less 
DO 10 1=1, NCLUHPS- 1 
DO 20 J*1,NCLUKPS-I 

IF (HIDTH(J).6T.NIDTH(J+1)) THEN 
CALL SNAP(NIDTH(J),NIDTH(J+U) 

CALL SHAP ( START ( J ) , START ( J+ 1 > ) 

CALL SNAP(SUH(J),SUH(J+1)) 

ENDIF 

20 CONTINUE 
10 CONTINUE 
RETURN 
END 
C 

SUBROUTINE SHAP(I,J) 

INTEGER I , J , TEHP 
TEHP=I 
I=J 

J=TEHP 
RETURN 
END 
C 

SUBROUT I NE COUNTCLUHPS ( CHAN , NCLUHPS ) 

INTE6ER HAXCHAN 
PARAHETER(HAXCHAN=140 ) 

INTE6ER CHAN , NCLUHPS , I , START ,HI DTH , SUN 
DIMENSION CHAN(0:HAXCHAN-1) 

NCLUHPS=0 
1=0 

10 CALL CLUHPINFOI CHAN , I , START , HI DTH . SUN ) 

IF (HIDTH.EQ.O) RETURN 
NCLUHPS=NCLUHPS+1 
I=START+NIDTH 

IF (I.LT.HAICHAN-1) 60 TO 10 
RETURN 
END 
C 

INTE6ER FUNCTION DUAL(N) 

INTE6ER N 
IF (N.LT.70) THEN 
DUAL=H0D(41IN, 701+70 
ELSE 

DUAL=HOD ( 4 1 IN , 70 ) 

ENDIF 
RETURN 
END 
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SUBROUTINE ZAPEND(NCLUHPS, CHAN, START, NIDTH, SUH) 

INTE6ER HAXCHAM 
PARAMETER ( HAICHAN=140 ) 

INTEGER I ,NCLUHPS,ENDPT, CHAN, START, NIDTH, SUN, DUAL, SSTART,NSUB 
INTEGER PREDICT 

DIMENSION CHAN(0:HAXCHAN-1),NIDTH(10),START(10|,SUN(10) 

DO 25 1=1 , HCLUHPS 
IF ( (HIDTHI I ) .6T.2) .AND. 

♦ (CHAN(START( I ) ) . GT . CHAM( START ( I ) *1 ) ) .AND. 

♦ ( CHANf DUAL ( START ( I ) ) ) . GT . 0 ) ) THEN 

C Extract the Gaussian shaped duap froa a possibly large duap 

C looting froa the left side of the duap 

CALL 6ETSUBCLUHP ( CHAN, START ( I ) , HIDTHI I ) , 6START , HSUB , 1 ) 

C Extrapolate the correct value for CHANISTART(D) 

CALL EXTRAPOLATE (CHAN , NSUB ,ST ART ( I ) *1 .START ( I ) .PREOICT ) 

C Adjust the aaplitudes of the endpoint and its dual 

sum i )=sum d-chanistarti i i ) ♦predict 

CHANISTART ( I ) )=PREDICT 

CHAN ( DUAL ( START ( I ) ) )=CHAN (DUAL (START! I) ) J-PREDICT 
C Readjust uidth and start pointers if necessary 

IF ( PREDICT. EB.O) THEN 
HIDTH(I)=HIDTH(I)*1 
START ( I )=START ( I )+l 
ENDIF 
END IF 

ENDPT=START(I)+NIDTH(I)-1 
IF ( (HIDTH( I ) .6T . 2) .AND . 

♦ ( CHAN ( ENDPT ) .GT. CHAN (ENDPT-1 ) ) .AND. 

♦ ( CHAN ( DUAL ( ENDPT U.6T.0) ) THEN 

C Extract the 6aussian shaped duap froa a possibly large duap 

C looting froa the right side of the duap 

CALL 6ETSUBCLUNP( CHAN, START( I ) ,HIDTH(I ) ,SSTART,NSUB,-1 ) 

C Extrapolate the correct value for CHAN(ENDPT) 

CALL EXTRAPOLATE (CHAN, HSUB, GSTART, ENDPT .PREDICT) 

C Adjust the aaplitudes of the endpoint and its dual 

SUH ( I ) =SUH ( I ) -CHAN ( ENDPT ) +PRED ICT 
CHANtENDPT )=PREDICT 

CHAN ( DUAL ( ENDPT )) -CHAN ( DUAL ( ENDPT ) ) -PRED I CT 
C Readjust uidth pointer if necessary 

IF ( PREDICT. E0.0) NIDTH(I)=NIDTH(I)-1 
ENDIF 
25 CONTINUE 
RETURN 
END 
C 

SUBROUTINE 6ETSUBCLUHPICHAN, START, NIDTH, 6START, HSUB, INCR) 

INTE6ER MAXCHAN 
PARAHETER(HAXCHAN=140) 

INTE6ER CHAN. START. NIDTH, 6START, HSUB, INCR.POS 
DIMENSION CHAN(0:MAXCHAN-1) 

C Find the 6aussian subduap starting froa the left if INCR=1 or 
C froa the right if INCR=-1 
HSUB-1 

C Set starting position for Gaussian duap search; endpoint is bad 
IF (INCR.E8.1) THEN 
P0S=START+1 
ELSE 

P0S=START+HIDTH-2 

ENDIF 
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C Find other end of 6aussian subciuap by diabing one hill until 
C the next hill is detected or end of given duap 
C First uphill 

10 IF (CHANtPOS) .LE.CHAN(POS«INCR) ) THEN 
NSUB»NSUB+1 
POSsPOSHNCR 

IF { (P0S.6T. STARTS) .AND. (P0S.LT.START+NIDTH-2) .AND. 

♦ (CHAN(POS) .61.0) ) GO TO 10 
EXDIF 

C Then domhill 

20 IF ( (P0S.6T. START) .AND. (POS.LT.START+NIDTH-1) .AND. 

* ( CHAN ( POS ) . 6E . CHAN ( P0S+ I NCR ) ) THEN 
NSUB*NSUB+1 

P0S«P0S*INCR 
60 TO 20 
ENDIF 

C Set starting pointer for the 6aussian sub duap 

IF (INCR.E0.1) THEN 
6START=START 
ELSE 

6START »ST ART *Nl DTH-NSUB- 1 
ENDIF 
RETURN 
END 
C 

SUBROUTINE EXTRAPOLATE ( CHAN ,NIDTH, XSTART , X SEARCH , PREDICT ) 
INTE6ER HAXCHAN , 0PP1 , 0PP2 , Y1 , Y2 
REAL C,X0PP, CENTER 
PARAHETER ( NAXCHAN=140 ) 

I NTE6ER CHAN , HI DTH .XSTART , XSEARCH, PREDICT 
DIMENSION CHAN(0:HAXCHAN-1) 

C Find the center C of the distribution 
C-CENTER ( CHAN » H I DTH , XSTART ) 

C Coapute the two x values closest to the opposite froa XSEARCH 
C These foraulas are correct ahether XSEARCH is to the left or the 
C right of C. 

0PPl s 2tC— XSEARCH 
0PP2=0PP1+1 

C 6et the opposite value froa XSEARCH relative to C as a real nuaber 
X0PP=2»C-XSEARCH 

C Set PREDICT to the linear interpolated value 

IF ( (0PP1.LT. XSTART). OR. (0PP1.6E.XSTART+HIDTH) ) THEN 
Y1=0 
ELSE 

Yl-CHAN(OPPl) 

ENDIF 

IF ( ( 0PP2.6E.X5TART+HIDTH). OR. (0PP2.LT. XSTART) ) THEN 
Y2=0 
ELSE 

Y2=CHAN(0PP2) 

ENDIF 

PREDICT= ( i . 01 Y2-Y1 ) / ( 1 . 0I0PP2-0PP1 ) » ( X0PP-0PP1 I +Y1 

RETURN 

END 
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REAL FUNCTION CENTER (CHAN, HI DTK, START) 

INTE6ER HAXCHAN, I, HIDTH, START, CHAN 
REAL YSUH,XDOTY 
PARAMETER ( HAXCHAN® 1 40 ) 

DIMENSION CHAN( Ot HAXCHAN- 1) 

C Deteraine the center of the duap starting at START of size HIDTH 
YSUH=0 
XD0TY=0 

DO 10 I=START,START+H1DTH-1 
YSUH=YSUH*CHAN(I) 

XDOTY=XDDTY+I ICHAN( I ) 

10 CONTINUE 

IF (YSUH.EQ.O) THEN 
NRITEIt, ERROR - YSUH = 0' 

CENTER»9999 

RETURN 

ENDIF 

CENTERsXDOTY/YSUH 

RETURN 

END 
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